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We develop a device model for p-i-n tunneling transit-time diodes based on single- and multiple 
graphene layer structures operating at the reverse bias voltages. The model of the graphene tunneling 
transit-time diode (GTUNNETT) accounts for the features of the interband tunneling generation 
of electrons and holes and their ballistic transport in the device i-section, as well as the effect 
of the self-consistent electric field associated with the charges of propagating electrons and holes. 
Using the developed model, we calculate the dc current-voltage characteristics and the small-signal 
ac frequency-dependent admittance as functions of the GTUNNETT structural parameters , in 
particular, the number of graphene layers and the dielectric constant of the surrounding media. It is 
shown that the admittance real part can be negative in a certain frequency range. As revealed, if the 
i-section somewhat shorter than one micrometer, this range corresponds to the terahertz frequencies. 
Due to the effect of the self-consistent electric field, the behavior of the GTUNNETT admittance in 
the range of its negativity of its real part is rather sensitive to the relation between the number of 
graphene layers and dielectric constant. The obtained results demonstrate that GTUNNETTs with 
optimized structure can be used in efficient terahertz oscillators. 



I. INTRODUCTION 



Pioneering papers by Shur and Eastman [T|, Q have 
stimulated extensive studies (which continue already for 
the fourth decade) of ballistic electron and hole trans- 
port (BET and BHT, respectively), i.e., collision free 
transport in short semiconductor structures. The main 
incentive is the realization of fastest velocities of elec- 
trons/holes and, hence, achievement of the operation of 
diodes and transistors in terahertz range (THz) of fre- 
quencies and low power consumption. Even at the initial 
stage of the ballistic transport research, several concepts 
of ballistic THz sources have been put forward (see, for 
instance, an early review Q). However, the realization 
THz generation in different semiconductor devices asso- 
ciated with BET/BHT, in particular, analogous to vac- 
uum devices, meets the problems associated with elec- 
tron scattering in real semiconductor structures. Cre- 
ation of heterostructures with selective doping with a 
two-dimensional electron gas (2DEG) spatially separated 
from the donors, has resulted in achievement of very long 
mean free path of electrons, at least at low temperatures. 
Recent discoveries of unique properties of graphene [3, [j| , 
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in particular, the demonstration of possibility of very 
long electron and hole mean free path in graphene lay- 
ers (GLs) and what is even more interesting in multiple 
graphene layers (MGLs) [(| 0| add optimism in build- 
ing graphene based THz devices using BET. The con- 
cept of graphene tunneling transit-time (GTUNNETT) 
p-i-n diode , which exhibits a negative dynamic conduc- 
tivity in the THz range, was proposed and substantiated 
in Refs. Q. This concept based not only on BET or 
quasi-BET (as well as BHT or quasi-BHT) in GLs and 
MGLs, but also on a strong interband tunneling under 
the electric field with a pronounced anisotropy [lCt ITl| 
due to the gapless energy spectrum, and constant ab- 
solute value of electrons and holes velocities Due 
to this, the electrons in the conduction band and the 
holes in the valence band generated owing to the inter- 
band tunneling in the electric field propogate primarily 
in the electric field direction with the velocity (in this 
direction) virtually equal to the characteristic velocity 
vw — 10 s cm/s. A large value of the directed veloc- 
ity in GLs and MGLs promotes the device operation at 
elevated frequencies. 

As shown [H, for the self-excitation of THz oscil- 
lations in a circuit with GTUNNETT diode, this circuit 
should serve as a resonator. However, at elevated tunnel- 
ing currents in GTUNNETTs considered previously |1] 
the self-consistent charge associated with propagating 
electron and hole streams can affect the spatial distri- 
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FIG. 1: Schematic views of GTUNNETT p-i-n diodes (a) 
with a single GL, (b) with an MGL structure, and (s) their 
band diagram at reverse bias. Arrows show the propagation 
directions of electrons and holes generated due to interband 
tunneling (mainly in those regions, where the electric field is 
relatively strong). 



bution of the self-consistent electric field and the electric 
potential in the i-section. As a result, the self-consistent 
electric field near the p-i- and i-n-juctions can be substan- 
tially reinforced. This, in turn, influences the tunneling 
generation of electrons and holes and the their transit 
conditions and, hence, the GTUNNETT dc and ac char- 
acteristics. 

In this paper, in contrast to the previous treat- 
ment 0, IH, we account for the self-consistent electric 
field associated with the variations of the electron and 
hole lateral charges in the i-section and their effect on 
the injection and the dc and ac characteristics. The 
effects of the space charge in planar TUNNETTs with 
the propagation of carriers perpendicular to the struc- 
ture plane were considered by Gribnikov et al. [l2j . The 
problems of calculation of the dc and ac characteristics 
of devices based on lateral structures accounting for the 
self-consistent electric field are substantially complicated 
by the features of the structure geometry (2D electron 
and hole channels and blade-like contact regions). In 
particular, as shown in the following, the related mathe- 
matical problems are reduced to a system of rather com- 
plex nonlinear integral-differential equations. Using the 
GTUNNETT device model, we derive these equations, 
solve them numerically, and find the characteristics. 



II. EQUATIONS OF THE MODEL 

The device under consideration is based on a GL or 
MGL structures with p- and n- side sections with the 
contacts and an undoped transit i-section. For dcfinitc- 
ness, we assume that p-and n-sections in a GL or in an 
MGL are created by chemical doping [l3j], although in 
similar devices but with extra gates these sections can 
be formed electrically [HI Hi} - Under the reverse bias 
voltage Vo, the potential drops primarily across the i- 
section A schematic view of the structures in question 
and their band-diagram under sufficiently strong reverse 
bias (Vo > 0) corresponding to the potential distribution 
in the i-section with a markedly nonuniform electric field 
are shown in Fig. 1. 

The electrons and holes injected into the i-section of 
one GL are characterized by the electron and hole sheet 
concentrations, £~ = £~(t,x) and S + = T, + (t,x), re- 
spectively and the potential (p = ip(t, x, y). Here the axis 
x is directed along the GL (or MGL) plane, i.e., in the 
direction of the current, the axis y is directed perpendic- 
ular to this plane, and t is the time. In the general case 
when the i-section is based on MGL, £~ = £~(t, x) and 
S + = S + (t, x) are the electron and hole densities in each 
GLs. In the case of BET and BHT, on which we focus 
mainly, the electron and hole sheet densities in each GL 
obey the continuity equations 
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respectively. Here e — |e| is the electron charge, 88 is 
the dielectric constant of the media surrounding GL (or 
MGL structure), K is the number of GLs in the GTUN- 
NETT structure, and S(y) is the delta function reflect- 
ing a narrowiness of GL and MGL structures even with 
rather large number of GLs in the y— direction. Equa- 
tion (1) corresponds to the situation when the electrons 
and holes generated due to the interband tunneling in 
the i-section obtain the directed velocities v x — vw and 
v x = —vw and preserve them during the propagation. 
The boundary conditions correspond to the assumption 
that the electrons and holes appear in the i-section only 
due to the interband tunneling (the injection of electrons 
from the p-section and holes from the n-section is neg- 
ligible) and that the highly conducting side contacts to 
the p- and n-sections are of blade type (the thicknesses of 
GL and MGL and the contacts to them are much smaller 
than the spacing between the contacts): 
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$|«<-i,v=o = -V/2, *|x>+j,»=o = V/2, (4) 

where 2/ is the length of the i-section and V = Vq + 
SVaj exp(-iujt) is the net voltage, which comprises the 
bias voltage Vq and the signal component with the ampli- 
tude SV and the frequency to. The inter band tunneling 
generation rate of electrons and holes in each GL (per 
unit of its area) is given by [§4Tl1| 



f2l\dip/dx\\ 3/2 
where <p{x) = $(x, y)\y=o, h is the Planck constant, 
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so that the characteristic tunneling dc current (per unit 
length in the transverse z-direction) and the characteris- 
tic electron and hole sheet density can be presented as 



Jni 



ielgo, 



So = 2lg /v w . 



From equation (2) with boundary condition (4), for the 
potential tp(x) in the GL (or MGL structure) plane we 
obtain 
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Introducing the dimensionless quantities: n T = 
E T /£ , £ = x /i } £' = x '/i 7 an d T = Vw t/l, Eqs. (1) 
and (6) can be reduced to the following system of non- 
linear integro-differential equations: 
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The term with the factor 7 in the right-hand side of 
Eq. (8) is associated with the contributions of the elec- 
tron and hole charges to the self-consistent electric field 
in the i-section. 

The dc and ac components of the terminal dc and ac 
currents, Jo and SJ^, are expressed via the dc and ac 
components and <5nJ as follows: 



J = ev w KE (n£ +n ), 



(11) 



6J u = ev w KE P (£)(dn+ + 6nZ)d£-iojC6V u . (12) 



Here p(£) = l/Vyl — £ 2 is the form factor and C ~ 
se/2-K 2 is the geometrical capacitance 001. The ex- 
plicit coordinate dependence of the form factor is a con- 
sequence of the Shocley-Ramo theorem 0, [Tj} for the 
device geometry under consideration. 

Hence the GTUNNETT small- signal admittance = 
SJuj/SVu is presented in the form 



ev w KE j p($—(6n++6nZ)d£-iojC. (13) 



As follows from Eq. (8), the problem under considera- 
tion is characterized by the parameter 7. 

If ae = 1.0, K = 1,2/ = 0.7 /im, in the voltage range 
V = 100 - 200 mV one obtains S ^ (0.6 - 1.7) x 
10 10 cm~ 2 and 7 ~ 0.61 - 0.87. 



III. SPATIAL POTENTIAL DISTRIBUTIONS 
AND CURRENT- VOLTAGE CHARACTERISTICS 

If the charges created by the propagating electrons and 
holes in the i-section are insignificant, that corresponds 
to 7 -C 1. In this case, the potential distribution is given 
by 
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and Eqs. (8), neglecting the term with 7, can be 
solved analytically. Taking into account boundary con- 
ditions Eq. (3), from Eqs. (8) we obtain 
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After that, using Eq. (11) and considering Eq. (15), one 
can find the following formula for the dc current-voltage 
characteristic: 



Jn = 



KJ, 
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FIG. 2: Spatial distributions of electron and hole concentra- 
tions (upper panel) and of electric potential (lower panel) in 
the GTUNNETT i-section. 



FIG. 4: Dependences of dc current on dielectric constant for 
different numbers of GLs K: net current - upper panel and 
current per one GL - lower panel. 



10,0 




50 100 150 200 

DC Voltage, mV 

FIG. 3: Current-voltage characteristics for GTUNNETTs 
with different number of GLs K. 
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FIG. 5: Dependences of dc current on number of GLs K at 
different voltages. 
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where T(x) is the Gamma- function. A distinction be- 
tween Jo and Joo is due to the nonuniformity of the elec- 
tric field in the i-section associated with the feature of 
the device geometry taken into account calculating Jo. 
At K = 1, 21 = 0.7 ^m and V = 100 - 200 mV, Eq. (16) 
yields J ~ 0.18 - 0.51 A/cm. 

To take the effect of the self-consistent electric field on 
the tunneling, one needs to solve system of Eqs. (8). Due 



to a complexity of the nonlinear integro-differential equa- 
tions in question, a numerical approach is indispensable. 
Equations (8) were solved numerically using successive 
approximation method, which is valid when 7 < j c — 5. 
In the cases 7 > 5, the method of the parameter evolu- 
tion was implemented. 

Using Eqs. (8) and (11) and setting V — Vo, we cal- 
culated numerically the GTUNNETT dc characteristics: 
spatial distributions of the dc electric potential and the 
dc components of the electron and hole densities, as well 
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as the dc current-voltage characteristics. The pertinent 
results are shown in Figs. 2-5. 

Upper panel in Fig. 2 shows examples of the spatial dis- 
tributions of the electron and hole sheet concentrations 
in the i-section. The spatial distributions of the dc elec- 
tric potential calculated for GTUNNETTs with different 
numbers of GLs K at fixed voltage are shown in Fig. 2 
(lower) panel. One can see that an increase in K results 
in a marked concentration of the electric field near the 
doped sections. This is because at larger K, the net tun- 
neling generation rate becomes stronger. This, in turn, 
results in higher charges of propagating electron and hole 
components, particularly, near the p-i- and i-n-junctions, 
respectively. 

Figure 3 demonstrates a difference in the dc current- 
voltage characteristics in GTUNNETTs with different 
numbers of GLs K. An increase in K leads to an in- 
crease in the dc current not only because of the increase 
in the number of current channels but because of the re- 
inforcement of the self-consistent electric field near the 
edges of the i-section(as mentioned above). As can be 
seen in Fig. 4 (upper panel), the three- fold increase in K 
leads to more than three-fold increase in the dc current 
particularly at relatively small dielectric constants. This 
is confirmed by plots in Fig. 4 (lower panel) , from which it 
follows that the dc current in each GL is larger in the de- 
vices with larger number of GLs. It is worth noting that 
when the dielectric constant is sufficiently large (about 
15-20), the dc current in one GL is virtually the same in 
the devices with different number of GLs, because in this 
case the charge effect is suppressed. In this case the value 
of the dc current (per one GL) is approximately the same 
as that obtained in the above numerical estimate using 
analytical formula given by Eq. (16). In GTUNNETTs 
with a moderate value of ae, the charge effect leads to 
supcrlincar dependences of the net dc current vs number 
of GLs as seen in Fig. 5. 

IV. GTUNNETT ADMITTANCE 

Disregarding the term in Eq. (8) associated with the 
electron and hole charges in the i-section, and using 
Eq. (12) , one can obtain analytically the following for- 
mulas for the ac potential distribution in the i-section, 
and the imaginary and real parts of the GTUNNETT 
admittance [9J : 

8rtx) u = ^ e w- 1 (j\ (17) 

lmY u = Y sin(wTi) • J {ur t ) - uC, (18) 

ReY u = Y Q cos(wr t ) • J Q {ur t ), (19) 

where Y = dJ /dVo — 3Jo/2Vb oc \/Vq is the dc dif- 
ferential conductivity, C is the geometrical capacitance, 
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FIG. 6: Frequency dependences of imaginary part, Im Y u , 
and real part, Re Y u , of GTUNNETT admittance (upper and 
lower panels, respectively) at different values of dielectric con- 
stant as (K = 1 and 21 = 0.7 fj,m). 
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FIG. 7: Frequency dependences of real part, Re Y u , 
of GTUNNETT admittance in the frequency range, where 
Re Yu < at different values of dielectric constant ae (K — 1 
and 21 — 0.7 /j,m) calculated accounting for electron and hole 
charges (solid lines) and neglecting them (dashed line) 



r t = l/vw is the characteristic transit time of electrons 
and holes across the i-section, and Jq{u)t) is the Bessel 
function. As follows from Eqs. (18) and (19), the admit- 
tance imaginary and real parts oscillate as functions of 
the transit angle urt with Re Y u < at ojTt near the tran- 
sit time resonances WTt = (2n—l/2)7r, where n = 1,2,3... 
is the resonance index. 
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FIG. 8: Frequency dependences of real part of admittance 
Re Yu, for GTUNNETTs with different number of GLs K (ae = 
10.0 and 21 = 0.7 /im). 



-5 

E 
E 

55 

E -10 



E. -15 



-20 



V„ = 200 mV 



■ 33 = 4.0 

• ae = 8.0 

as = 12.0 



1 



5 10 15 20 25 30 35 40 
Number of GLs K 



time resonances. As seen from Fig. 7 , the minima of 
the admittance real part (where the latter is negative) 
become deeper and the minima shift toward smaller fre- 
quencies with increasing dielectric constant ae. This is 
because an increase in ae leads to weakening of the ef- 
fect of the electron and hole charges on the injection and 
propagation processes. The dashed line corresponds to 
the case when the role of these charges is diminished (very 
large ae) considered previously @ . 

Figure 8 shows the admittance real part vs signal fre- 
quency calculated for GTUNNETTs with different num- 
ber of GLs K. One can see that GTUNNETTs with 
larger number of GLs K exhibit much more deep min- 
ima (compare the curves which correspond to K = 1 and 
ai = 10.0 in Fig. 7 and that for K = 10 and ae = 10.0 
in Fig. 8.). As follows from this comparison, the ten- 
fold increase in K leads to about four-fold increase in 
the minimum depth. 

Figure 9 shows the dependence of the minimum value 
of the real part of the GTUNNETT admittance min 
{Rel^j} on the number of GLs K calculated for differ- 
ent dielectric constants ae. The nonmonotonic character 
of these dependences is associated with the interplay of 
the factors determined by the influence of electron and 
hole charges: an increase in ae leads to a suppression 
of these charges role, while an increase in K results in 
proportional reinforcement of the effect of these charges 
accompanied by proportional increase in the current. 



FIG. 9: Minimum value of real part of GTUNNETT admit- 
tance as a function of number of GLs K for GTUNNETTs 
with different dielectric constant ae. 



The numerical solution of Eqs. (8), accounting for the 
ac components of the applied voltage 5V U and the the 
self-consistent charges and electric field, shows that the 
oscillatory frequency dependences are preserved in this 
more realistic case although they are be quantitatively 
modified. In particular, these dependences vary with 
varying such parameters as the number of GLs K and 
dielectric constant ae. The results of numerical calcula- 
tions using Eq. (8) are demonstrated in Figs. 6-9. 

Figure 6 shows the dependence of the imaginary and 
real parts, Im Y^ and Re Y m of the admittance on 
the signal frequency calculated for GTUNNETTs with 
21 = 0.7 /im and K = 1 at V = 200 mV at different val- 
ues of dielectric constant ae. As follows from Fig. 6, the 
imaginary and real parts of the GTUNNETT admittance 
are oscillatory functions of the signal frequency [see also 
Eqs. (18) and (19)]. The oscillatory behavior is due to 
the transit-time resonances. At sufficiently high frequen- 
cies, the admittance imaginary part is negative because it 
is determined primarily by the geometrical capacitance, 
so that Im Y u ~ — uj C. In the certain frequency ranges, 
the admittance real part is negative as shown in Fig. 6 
(lower panel). These ranges correspond to the transit- 



V. DISCUSSION 

Calculating the dc current, we disregarded the contri- 
butions of the thermogeneration of electron-hole pairs in 
the i-section and the injection of minority carriers from 
the n- and p-sections. 

The thermogeneration dc current is given by 

jtherm = 4 Kdg therm_ ^ 

The thermogeneration rate g* ) herm at the temperature 
T = 300 K is estimated as [H gf er ~ 10 20 - 
10 21 cm _2 s -1 . As a result, at K = 1 and 21 
0.7 fxm from Eq. (20)one obtains J* herm ~ 2 x (10^ 3 - 
10~ 2 ) A/cm. 

The injection dc current can be estimated as 

inj _2KeT^ (_t\ (21) 

where ji is the Fermi energy of electrons in the n-section 
and holes in the p-section. Setting ji = 50 meV, from 
Eq. (21) for K = 1 we find ~ 0.077 A/cm. Thus, 
jtherm^ j™3 ^ & ^ least in the voltage range consid- 
ered above. 

The typical values of the real part of the GTUNNETT 
small-signal admittance in the vicinity of the first transit- 
time resonance (wr t ~ 3/2) is in the range Re Y u ~ —5 to 
—20 mS/mm depending on K and as. This implies that 
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the ac resistance, R u = H/Y omega of optimized GTUN- 
NETTs with Re Y u = -20 mS/mm (K = 10, ac = 10, 
and V = 200 mV) and the width H = (0.5 - 1) mm, 
is R u ~ 50 — 100 SI. Assuming that in this case 
J ~ 0.6 - 0.8 A/mm at Vb = 200 mV, the dc power, 
Po, generated by the propagating electrons and holes in 
the contact n- and p-regions of the device can be esti- 
mated as Po — 150 mW. When Re Y u = —5 mS/mm 
(K = 1 and ae = 10), the same resistance is provided if 
H = 2 - 4 mm with the dc power P ~ 20 - 40 mW. 

The depth of the minimum of the real part of the 
GTUNNETT admittance can be increased by applying 
higher bias voltage. This is because the admittance is 
proportional to the dc differential conductivity which, in 
turn, is approximately proportional to y/Vo- However, 
an increase in Vq results in the pertinent rise of Pq. If 
eVb exceeds the the energy of optical phonons in GLs 
Huj ~ 200 meV, their fast emission disrupts BET and 
BHT in the i-section. The emission of optical phonon 
leads to some isotropization of the angular distribution 
of the electron and hole velocities of electrons and holes. 
This results in some decrease in their mean directed ve- 
locities (they become somewhat smaller than ivw) and, 
hence, in some increase in the electron and hole charges 
and the transit time. Such factors promote more steep dc 
current-voltage characteristics and modify the frequency 
dependence of the admittance. Although they should not 
change the GTUNNETT performance qualitatively. 

The depth of the first resonance minimum of the real 
part of the GTUNNETT admittance with the param- 
eters used in the above estimate is of the same order 
of magnitude as that of the small area (about 2 /itm 2 ) 
THz resonant-tunneling diodes considered by Asada et 
al. [l9j . It is worth noting that the frequency at which 
Im Yu turns zero can fall into the frequency range, where 
Re Y u < 0. In such a case, GTUNNETTs can exhibit 
the self-excitation of THz oscillations even without the 
external resonant circuit. Moreover, the contact n- and 
p-sections (both ungated as in the case of their chemical 
doping or gated in the case of the electrical "doping") 
can serve as plasma resonant cavities [9|, [15[ . The com- 
bination of the transi-time and the plasma resonances, 
can substantially liberalize the self-excitation conditions 
and enhance the THz emission efficiency. The match- 
ing of the transit-time and plasma resonant frequencies 
requires a proper choice of the length of the n- and p- 
sections (from several micrometers in the ungated struc- 
tures to about one micrometer in the gated ones. The 



plasmons and related resonances in the gated GLs were 
predicted and analyzed several years ago [2(],[2l| (see also 
Ref. 22] . Recently, the plasmons in gated GLs were de- 
tected experimentally [2J, [24| . It is interesting, that as 
a gate for the GL plasmon cavity another gate can be 
used [111 when the contact regions constitute double-GL 
structures (26|.[27|. 

VI. CONCLUSIONS 



We have developed a self-consistent device model for 
p-i-n GTUNNETTs with different number of GLs which 
enables the calculation of their realistic dc and ac char- 
acteristics. Our calculations have shown that: 

(i) The charges of electrons and holes propagating in the 
i-section and their effect on the spatial distribution of the 
self-consistent electric field increase the steepness of the 
dc current-voltage characteristics. This effect is stronger 
in GTUNNETTs with larger number of GLs and can be 
weakened by using the media surrounding GL or MGL 
structure (substrates and top dielectric layers) with ele- 
vated values of the dielectric constant; 

(ii) The imaginary and real parts of the GTUNNETT ad- 
mittance are oscillatory functions of the signal frequency 
due to the transit-time resonances. In the certain fre- 
quency ranges, which correspond to the THz frequencies, 
the admittance real part can be negative that enables the 
use of GTUNNETTs in the sources of THz radiation; 

(iii) The charges of electrons and holes influence on the 
spatial distribution of the self-consistent electric field and 
the interband tunneling and can substantially affect the 
GTUNNETT admittance, particularly, its frequency de- 
pendence near the transit-time resonances. The role of 
this effect enhances with increasing number of GLs in 
the device structure and can be markedly weakened in 
the devices with relatively large dielectric constant; 

(iv) The sensitivity of the GTUNNETT characteristics 
to the structural parameters (number of GLs, dielec- 
tric constant, length of the i-section) and the bias volt- 
age opens up wide opportunity for the optimization of 
GTUNNETTs for THz oscillators. 
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